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Abstract 

We present the two-loop virtual QCD corrections to the production of heavy quarks in the quark- 
anti-quark-annihilation channel in the limit when all kinematical invariants are large compared 
to the mass of the heavy quark. Our result is exact up to terms suppressed by powers of the 
heavy-quark mass. The derivation is based on a simple relation between massless and massive 
scattering amplitudes in gauge theories proposed recently by two of the authors as well as a direct 
calculation of the massive amplitude at two loops. The results presented here form an important 
part of the next-to-next-to-leading order QCD contributions to heavy-quark production in hadron- 
hadron collisions. 



The hadro-production of heavy quarks, especially of top-quarks, is an important process at 
hadron colliders. Thus far, the Tevatron has provided us with a wealth of information on the top- 
quark, most prominently with a precise measurement of its mass. In the future, the LHC is expected 
to accumulate very high statistics for the production of tt-pairs, approximately 8-10^ events per 
year in the initial low luminosity run [1]. At the LHC the uses of ?f-pairs are e.g. for energy 
scale calibration, for background estimates and last but not least, for precision measurements of 
Standard Model parameters. 

The present and in particular the anticipated future experimental precision on heavy-quark 
hadro-production requires theory predictions to include radiative corrections in Quantum Chromo- 
dynamics (QCD) beyond the next-to-leading order (NLO). For instance, at the LHC the total cross 
section for ^/"-production at NLO in QCD is accurate to 0(15%) only, which has to be contrasted 
with expected precision for the top-mass measurement of o(lGeV). Another important process 
is the production of Z?-quarks at moderate to large transverse momentum. NLO QCD corrections 
for heavy-quark hadro-production have been known since long, see e.g. Refs. [2-8]. Knowledge 
of the radiative QCD corrections at next-to-next-to-leading order (NNLO) will certainly improve 
the stability of theory predictions with respect to scale variations and provide a match with precise 
parton evolution at NNLO [9, 10]. 

In this letter, we present results for the virtual QCD corrections at two loops for the pair- 
production of heavy quarks in the ^(^-annihilation channel. To be precise, we calculate the inter- 
ference of the two-loop amplitude with the Bom one and we work in the limit of fixed scattering 
angle and high energy, where all kinematic invariants are large compared to the heavy-quark mass. 
Thus, our result contains all logarithms in the heavy quark mass as well as all constant contribu- 
tions (i.e. the mass -independent terms). Throughout this letter we neglect power corrections in the 
heavy-quark mass. 

In our calculation we employ two different methods. On the one hand, we apply a general- 
ization of the infrared factorization formula for massless QCD amplitudes [11, 12] to the case of 
massive partons [13]. In essence, it results in an extremely simple universal multiplicative relation 
between a massive QCD amplitude in the small-mass limit and its massless version [13]. In this 
way, we can largely use for our derivation the results of Ref. [14], where the NNLO QCD correc- 
tions to massless quark-quark scattering (i.e. qq — > q'q') have been computed. On the other hand, 
we perform a direct calculation of the relevant Feynman diagrams in the massive case followed by 
a subsequent expansion in the small-mass limit. As an added benefit, this approach provides us 
with a non-trivial, albeit complicated, check of the massless results. Moreover, it makes it possible 
to systematically calculate power corrections in the mass, which can improve the convergence of 
the small-mass expansion. This is certainly relevant in the case of the top-quark pair-production at 
the LHC (less so, perhaps, for Z?-quark production). 

We would like to emphasize that the agreement between our prediction based on the approach 
of Ref. [13] and our direct calculation constitutes the first non-trivial check of this factorization ap- 
proach at two loops. The formalism of Ref. [13] has recently been also applied to the re-derivation 
of the two-loop QED corrections to Bhabha scattering [15] confirming earlier results on the two- 
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loop radiative photonic corrections [16-19] in that process. 



Setting the stage 

The pair-production of heavy quarks in the ^^-annihilation channel corresponds to the scattering 
process, 

q{Pl)+^{P2) Q{p3,m)+Q{p4,m), (1) 

where p, denote the quark momenta and m the mass of the heavy quark. Energy-momentum 
conservation implies 

^1+^2=^3+^4- (2) 

Following the notation of Ref. [14] we consider the scattering amplitude 5Vf for the process ([T]) at 
fixed values of the external parton momenta pi, thus p\ = P2 = ^ and pi = P4 = m^- It may be 
written as a series expansion in the strong coupling a^, 
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(3) 



where we define the expansion coefficients in powers of CXs(a'^) / (2;t) with n being the renormal- 
ization scale. We work in conventional dimensional regularization, J = 4 — 2£, in the MS-scheme 
for the coupling constant renormalization. The heavy mass m on the other hand is always taken to 
be the pole mass. 

We explicitly relate the bare (unrenormalized) coupling to the renormalized coupling tts by 



(4) 



where we put the factor 5e = (47r)^exp(— e'jt) = 1 for simplicity and (3 is the QCD |3-function 
[20,21] 



11 2 „ 17 2 5 

j2 



|3o = — C4--7>n/, ^i = —CA--CjJpnf-CpTpnf. (5) 



As we work in a general non-Abelian SU(A'^) -gauge theory we set = N, Cp = {N — \ )/2N 
and Tp — 1/2. Throughout this letter, denotes the number of colors and Hj- the total number of 
flavors, which is the sum of light and n^^ heavy quarks. 

The squared amplitude for the process ([T]) summed over spins and colors is a function of the 
Mandelstam variables s, t and u given by 

s = {p\+P2) , t = {pi-p3) -m , u = {p\-pa) -m . (6) 

Then it is convenient to define the function J?(£,m,5,?,/j) for the spin and color averaged ampli- 
tudes as 

Y,\^{q + q^Q + Q)? = -^Sl{z,m,s,t,^), (7) 
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which has a perturbative expansion similar to Eq. ([3]), 

Ji{£,m,s,t,/S) = 167i^a; 
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(8) 



In terms of the amplitudes the expansion coefficients in Eq. ([8]) may be expressed as 



t^ + u' 



£ + (m) , 



(5Vf(0)|^(0))=2(iv2_ J 

((fAf(0)|f^(l)) + (fAf(l)|fAf(0)) 
(<Af (l)|5Vf (1)) + (fAf W |5Vf (2)) + (2) (0)) 



(9) 
(10) 
(11) 



where we have neglected powers in the heavy-quark mass m in Jl^. Expressions for Sl^ with 
the complete heavy-quark mass dependence using dimensional regularization can be obtained e.g. 
from Ref. [7, 8]. The loop-by-loop contribution {M ^^^) in dimensional regularization in J^^ 
and also with the full heavy-quark mass dependence can be computed with the help of Ref. [22]. 
In this letter, we provide for the first time the real part of ^^^) up to powers (m) in the 

heavy-quark mass m. 



Massive amplitudes from QCD factorization 

Let us briefly recall the key features of Ref. [13] to calculate loop amplitudes with massive partons 
from massless ones. The QCD factorization approach rests on the fact that a massive amplitude 
£W [p]'('") for any given physical process shares essential properties in the small-mass limit with the 
corresponding massless amplitude !M [p]'('"=0). The latter one, [p1'('"=o)^ generally displays two 
types of singularities, soft and coUinear, related to the emission of gluons with vanishing energy 
and to coUinear parton radiation off massless hard partons, respectively. These appear explicitly 
as factorizing poles in £ in dimensional regularization after the usual ultraviolet renormalization is 
performed. In the former case, the soft singularities remain in 9\{ [p1'('") as single poles in £ while 
some of the coUinear singularities are now screened by the mass m of the heavy fields, which gives 
rise to a logarithmic dependence on m, see e.g. Ref. [23]. 

Thus, in the small-mass limit the differences between a massless and a massive amplitude can 
be thought of as a mere change in the regularization scheme. As an upshot, QCD factorization 
provides a direct relation between £W IpM"') and [p]'("'=o) which can be cast in the remarkably 
simple and suggestive relation 

^IviH = Yl (zj™l°))' X fAf[p]'('"=o). (12) 

I'e {all legs} 

The function Z^'"!^) is process independent and depends only on the external parton, i.e. quarks 
in the case at hand. For external massive quarks Q it is defined as the ratio of the on-shell heavy- 
quark form factor and the massless on-shell one, both being known [24-26] to sufficient orders in 
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tts and powers of £. An explicit expression for 

up to two loops is given0in Ref. [13]. Exploiting the full predictive power of the relation Eq. (fT2l) 
and applying it to the process Eq. ([T]) we get 

which assumes the hierarchy of scales <^ s,t,u , i.e. we neglect terms 0(m). Eq. (fT4)) predicts 
the complete real part of the squared amplitude (iW (^^) except for those terms, which are 
linear in n^^, i.e. the number of heavy quarks. These two-loop contributions have been excluded 
explicitly from the definition [13] of Z^'"!^^, as one needs additional process dependent terms for 
their description. Their incorporation for the case of Bhabha scattering was presented in Ref. [15] 
in agreement with the direct calculation [27]. 



Direct calculation of the massive amplitude 

An alternative method is the direct calculation of all necessary massive Feynman diagrams together 
with an expansion in the small mass. The advantage of this approach is an independent check of 
Eq. ([14)) as well as of the corresponding massless results. Moreover, it also allows for a relatively 
easy access to all heavy-quark loop corrections. Presently, these have not been obtained from the 
QCD factorization method, since they are related to the process dependent contributions. 

We performed our computation using the DiaGen/IdSolver system of one of the authors (M.C.). 
After the diagram generation phase, all the integrals have been reduced to a set of 145 masters with 
the help of the Laporta algorithm [28] extended by topological symmetry properties, where, how- 
ever, integrals which are related hy at ^ u exchange of the Mandelstam variables are considered 
as being independent. The evaluation of the masters proceeded similarly to the methodology de- 
veloped in Ref. [29,30]. The main idea was to construct Mellin-Bames [31,32] representations for 
all the integrals, followed by a subsequent analytic continuation with the MB package [33] and an 
expansion in the mass by closing contours. The resulting integrals have then been transformed into 
series representations, some two-fold, and resummed with the help of XSummer [34] in the cases 
of non-trivial dependence on the kinematic variables. Some constants, though, were not given by 
harmonic series and have been computed with the help of the PSLQ algorithm [35]. Needless to 
say that all of the above steps were performed fully automatically. 

In our calculation, we have only been able to obtain the leading color term and the full depen- 
dence on the number of light and heavy fermion species. The reason is that the terms subleading in 
color contain contributions from non-planar graphs. Fig. [T] shows the bottleneck cases with 6- and 
8-fold Mellin-Bames representations, for the six and seven liners respectively. The representations 

^Note the different normalization as/(47i) of the coupHng used in that reference. 
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Figure 1: The most complicated non-planar topolo- 
gies, together with the multiplicities of the respective 
master integrals for qq QQ scattering. The thick 
lines are massive. 



Figure 2: A heavy-quark loop 
contribution to the wave-function 
renormalization of the light quarks. 



have been obtained similarly as in Ref. [32], i.e. from the Feynman-parametric representation of 
the two-loop graphs and not by integrating loop-by-loop, the strategy adopted for planars. In this 
way, the mass expansion generates integrals that are at worst as complex as the massless ones, plus 
some new terms that behave as 1 / Vm^ and must cancel in the complete result. The sheer amount 
of remaining difficult 4-fold integrals turned out to be presently intractable with the developed 
software. 

Finally, as mentioned above, we renormalized the mass of the heavy quark and the external 
wave functions in the on-shell scheme (besides the MS renormalization of the coupling constant). 
The former results are now known up to three-loop level from Ref. [36, 37]. A peculiarity of the 
light quark states is that they obtain contributions to the wave-function renormalization constants 
from heavy-quark loops, see Fig.|2l In fact, we have 



Z2 



1 



4£ 2 



m 



5 ■ 
24 



(15) 



Results 

We are now in a position to present our result for qq QQ scattering for the interference of the 
two-loop and Bom amplitude, 

2Re(fW(°)|?if(2)) = (16) 
2{N^-\)(n^A + B + ^C + NniDi + Nn^Ph + + + (n, -f , 

which we choose to express by grouping terms according to the power of the number of colors 
and the numbers of light and = 1 heavy quarks with Hj- = + nj^ total flavors. As detailed 
above, the coefficients A, D/, Ei and F have been computed with both methods, i.e. by employing 
our universal multiplicative relation (fT2l) between the massive and massless amplitudes as well as 
a direct evaluation of the loop integrals in the small-mass expansion. The terms linear in rif^, 
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and Eh could also be easily obtained from a direct Feynman diagram calculation, while for B and 
C the approach based on factorization proved more powerful. 

As the only dimensionless kinematic variable in our problem, we choose x = —t/s and keep 
the dependence on the renormalization scale, /j, explicit. We also introduce the following compact 
notation 



log 



m 

s 



Lx = log(jc), L^ = log(l-j£:) . 



(17) 



The different components now read 



1 



[x^ X n If 


" 2 1" 


+u 


2 1' 


2lx^ 




X + - 


-x'+x--^ 


+ 4 











+2x - 1 ) + y I + 1 |l„,L, [-2x^ + 2x-l]+ 



1 



21a; 
29x2 



+ L,{-2x' 
29x 



+L^{-2x^ + 2x- l) + 



23' 
12 



19^2 I9x 



26x^ 55x 23 \ 173^2 
x^ X 

-y+3- 



+ L,(4x2-4x + 2) 
173x 



72 72 
+ L,„L^ [2x2 - 2x + l] + 



13 
12 

2 



+ (2x2 



+ ^ '-6+6-12 



2x2 2x 1 
'~ + y- 3 



-x^+x + lJx^-x + 



47x 35 



X X 



-4x2 



+■ 



12 8 



-4x^ + 5x - 



-4x-2) -- 
5 



3" 



7x2 1 
._ + _ + L,(4.^-4. + 2)-- 



47x2 



8x2 



Ux 1\ 487x2 
- + 3^^ + ^ 



487x 
~36~ 



,,x2 X 1\ 601 

+'■17-3 + 6+^ 



5x 



-^ + 2-4'+-'^ 



+Li2(x^ 

9907x2 
432 



1 



-2x ~ 2 1 ^ 



43x2 



9907x 2 



15 Ix 
25 



o , 4x2 5 .^. 5 



1 



23x2 



17x 17\„ 





"x2 X 1 ■ 




"2x2 


2x r 






7-4 + 8 






-y+3 


+ L 



4x2 



+- 



4x 





x2 X 1 






y-3 + 6 





11x2 llx 

l8~ + T8"+^-'^" 



2 

X X 

T + 3 



36 



+ y + L,. (-2x2 + 2x- 1) _ 5 
8x2 4 

~-T + 3 



+ L L2 



— + y + L,(-4x2 + 4x-2) 





" 14x2 / 




[9-9 





10 

"9" 



+ L2 



1 



3x 



247x2 247x 283 
+^6 36~ + ^ 



1 \ T 7 2 . \. 23x2 23x 83 
x--]lI+{2x^-3x)U + ^-- + - 



6 



35 



+ 



Ux 10\^ 



37x 37x 9 
—^ + -^ + 71 



+ L„ 



-4x2 + 



x^L^ 



+ 



19x 



■ + 



1 



2 4 
2 |L 



1 



4 

781^2 
72 



781x , 



Ix' 



lOx 



+Li3(x) (2x^-x+-] + {- ^ + tK3-T7T 



499 



144 



72 



7x2 



2 

X X 

T + 3 
1 



2x +x--]Ljc 



X 

12^ 3 



1 

24 



+ (^^ + L, (^2^2 - X + + l2 + Li2(x) (4x2 _ 2^ + 1) + 



^-2x+i)l,3 

86x2 j73^ 
— ^ + ^— 



+71^ 



+71^ 



8x2 5^ 



+ 



43x2 61x 91 \ 

36~^36"~72y ^ 

7x2 ^^0x2 

^ h Lv 

18 12 



49 

T J~-' ' 216 
23x2 



2003x2 2003x 



216 



+ Li3(x) (-4x2 + 2x- l) 



3 

2x+ 1 



17x 17\ ^ 919 



12 



13x 13\ lOlx . 
+1^+12 + ^ + " 



432 

x2 

' 2 
3 



+ 
5x 



X 1 , 
— + — L' 
24 48 



+ ■ 



11 



7x2 



25x2 25x 25^ 
"~6~'^1T~24 j 1-x 



^)l2 + S,,2(x)(-2x2 



11 , 

-5x+y ]U + 



260x2 139x 



+ 



52x2 32 
^ 3~^"6" 



^3 + 



4 

1027 \ 
72 J 



+ 71" 



25x2 53x 47 



+ — - — +L,.7l^ + ^ - -TT 



L, + . 



36 18 
17845x2 17845x 



llx 17 



2592 



2592 



+ 71-* 



11x2 

~80" 



12 
X 1 
240 "^480 



r-r^f^2 o 3\ ^ , , 2 . 11\ 2/^1009x2 259x 403\ 
+m(x) ^2x2 + 3x - - j + S2,(x) (2x2 + 5x - - j + 7^2 - ^ + _ j 

, , / 7x2 13x / 23\ 

+Li3(x) ( — ^ ^+U(-8x2 + 2x-l)+L,, (2x2 + 5x-y j - — j +Li2(x)((7x 



7x 7 

■y+4 



L? + 



7x2 23^ 



-2x2 _ 5;^- _^ . 



11 



13 \, 2 

+ — |Lv + 7l2 



■ + ■ 



17 
12 



+Lv 



-2^2 



5^ + y K3 + 



14x2 55x 149 \ 77287 
~9~ ^ Ts" ^ ^ 5184 



^ |-x2 +x- ^1 + ^ [-2x2 + 2x- l] [2x2 -2x+ l] - + (-4x' 

27 1 1 

+4x-2)+U(6x2-6x + 3) - y ^ + p-{L;„L, [4x^-4x + 2] +l2 [-2x2 + 2x- l] 



47x2 47x , . ^ T 7 . .\ 35 

- + + ( -4x2 + 4x - 2) + ^^^^2 _ + 3) _ _ 



+ L, 



49x2 49x 



+U (- 12x2 + 12x - 6) + Lv (8x2 _ 8^ _^ 4^ _^ 
143x 



29 \ 

+ Lj.(4x2-4x + 2) + y 



37 
12 

415x2 415x 



+ -6x2 ^ . 



15x 15 



Lf + 



67x2 



36 



■ + ■ 



36 



+ L, 



52x2 49^ 20' 



7 



+71" 



nx 17 

■ + ■ 



12 12 ' 24y ' >' 

3 



+ Lt 2x^-x + - 



17 
9" 



+ ■ 



2^2 2x 1 



+L,„l2 [-4^2 + 4x - 2] + L; 
+Ly (2x2 - 2x + 1) + 1] + L,„L, 



"4x2 


4x 2' 






~Y 


-T + 3 




2x2 _ 2;^; _|_ L -^ ^ 



25x2 25x 



+ U (- 12x2 + 12x - 6) + Ly (8x2 + 4) 



13 

+ — 



9x2 5 

-— + — + Lj,(-8x2 + 8x-4) +Lv(l2x2- 12x + 6) -- 



9x 



--)L^+(3x2-lf )L, 



7x 7 \ 5 
T + T2)+4 



2 

+ L, 



1 



16x 



/3x 



7x2 



x--j + — + L,.(-2x^+x+l)+7i^( — 



12x2 -15x + ^ 1^2 + 



46x2 55x / „ 9 „ ,N 
-^ + — +Lv(-8x2 + 8x-4) 



14\ 85x2 85x . , 2 ^ . 
-y U + -^-- + L-(-4x2 + 2x-l] 



17x2 j7 

^ + ^-T2l+^>' 



16x2 



lOx 

3~ 



31 
18 



+ 



15x 



3x+- — + L,, 3x 



X 1 
2 + 4 



+Li2 (x) 6x2 - 3x + - _^ 



2 ) 



50x2 137^ 

^ h 71 

3 12 



31x2 4^ 
^ + ~6~~ 12 



11 



2x2 



7x 7 



+- 

( 

2 
~3 
1 

~3 
1 

~6 



3 3 

47 
+ — 



+ ■ 



+ 71^ 



/ 263x2 lOlx 91 



V 72 



72 144 



+ L, 



86x2 j93^ 



+ 71" 



8x2 



+ 



91x2 73\ 4^3^ 

+ — I ^3 + ^ 



6 12y ^ 108 j 



+ L! 



2 

X X 

-2+2 



+ L,„L,^ 

+ L;, 



8x2 4' 
~~y + 3 



2x2 2x I 
'~ + y ~ 3 



+ L 



11x2 ^j^^ 



18 



llx 

4x2 4^ 

■- + T 

2x2 2x 



-x+- 



1 



1 

36 



+ l2l 



2 

- y + I + Lv (-4x^ + 4x - 2) + U (6x2 _ 6^ + 3) 



+ LmL^ 



I4x^ I4x T o 1 

+ ~ + (~^^ + 8x - 4) + U- (12x2 - + 6) - - 



+ L? 



16x 



-— + L4-8x" + 8x-4)+L,( ^ _ + _!+_ 



9 
3x2 
' 2 
229 
12 



16x2 



16x 8 



3x 



16x2 



9x 



157x2 ^ 2 

Lv hL: 

36 



4 2 



+ ■ 



157x 
36" 



+ 71" 



2x2 2x 1 
-~ + y ~ 3 



+ L, X 



+ L,j;L_5 



X _ 1 

2 ~ 2 

,2 



^ - 3x^ l2 + (9x - 6x2 ) ^ ^ l2 ( 1 - 2x) - 7x + 7i2 ^ - ^ + y 



-Z^+L,(4x2-2x-2)-^ 



+ L2 



-12x2 + 15x- 



15 



)l? + 



2x2 



llx 



+ L„(8x2-8x 



-4) 



, Lv + 
3/ 9 



40x2 40x 2/^17x2 17x 17 \ 



+ 71- 



+ ^j+Lv(4x2-2x+l)+L, 



28x2 34x 



26\ 161 

'T J +^ 



+ L„ 



-x2l^ + (I + L, (3x2 - I + ^) - ^) l2 + Li2 (x) (6x2 _ 3^ 



8 



+^)L,+ (m2_57^^^2 



-3a;2 + — 



6 U 



115x2 . 



+ ■ 



115x 



+U3{x)l-6x^ + 3x--\ + Siaix) {-4x'^ + 6x- 3) +lU —- — + -] +7^1 —-x 



24 



+L, 



+ L,.( -Sx^ + ^+Ti^{2x^-3x + 



+ 



37^2 



2Sx 14\ ^ 



67 

Is 



+ 6x - 3) + ^ + L„ (^-6^2 + X - + + Li2(x) ( - + 6x - 3) L, 



+ (^(2x-1)l2 + 2L,,+ 
2x\ 1957x 



100^2 181x 



+ ■ 



54 



3 

4x2 



+ 71^ 



62x2 4j 



3+-3 



14x 7 



+6) + Lij (x) ( 1 2x2 - 6x + 3) + ^ 

+ 



3 

172x2 



+ r 



/67x2 229x ^ 42l\ 
\36 36" 



1957x2 



54 



■ + ^(3 



+ +Si,2(x)(8x2-12x 



91x2 73 \ 553 ■ 



■ + 57x + 



1 



16x2 22x 



11 

y 



97 
y 



8 16 



35x 



18 12 



32x2 



+- 



16x 8 



7 



/ 7x2 5-^ 17 V T 



29x 13 

i2"~y 



365x 
72" 



95x2 



103x ^ 103 



12 



24 J 1 



7 1027 
~^ 144 



2x2 4^ 2 



L^+((-— +--Tm+(2-4)L;+(^ 



43x2 



+- 



15x 11 



2 L, 



304x2 J 53^ 

+ ■ 



1181 
^72~ / 



+ 



9 4 
38959x2 38959;,- 



125x2 149x 155 
+ ■ 



18 



36 ^ 36 



+ (-54x2+45x 



648 

25x 25 



648 



+ S2,2W (^-10x2 + 25x-^^ +Li4(x)(-6x2 



x + - +U -x-" + +L, + l+it 

o ; > I 12 24 n n ifi I 



7x2 



2x 13 



36 



+Si,3(x) (4x2 + 2x- _^l2 
/ 14x2 47^ 



28 Ix 



+ 71^ 



-4x2 



19x 



19 



77x2 29x 101 \ 
2?'^240y 

163 5 

~T2^x 



+S 



1,2 



+ L, (8x2 _ + 4) + L, (^10x2 - 9x - + y ^ + Li3 (x) 



+^ + L, 10x2 - 3x + y ^ + L, (24x2 _ + 5) + + Li2(x) (^-21x2 + ^ 



21 



23x 



15 



29 



+ 



43x2 28x 169 

yy^y is" 



,„ ^ /520x2 



2 

829x 



25 



6 12 



50x2 28x 25 

"y~^y~~y 



+ 



98x2 



+■ 



107x 71 \ ^ 605 \ y 
41473 



9 18 

3797x2 1585x ,^ . . ,^, — \ 

\ hf2x2-2x+l)log(2) 

216 108 ^ ^ J B\ J J 



3 

2353 \ 



1296 



(19) 



9 



c 



I jx'^ X 






" 2 1' 




2 1" 








. -x+- 




-X +X-- 













5x^ 5x 

H hLv 

2 2 



-4;c^ + 4x 



-2) + Lv (4;c^ -4x + 2) + 1} + |l„,L, [-2^;^ + 2x - l] + - ;c + ^ + L„, [3^;^ - 3;c 

+U (-4x2 + 4x - 2) + Lj, (4x2 _ 4;^. + 2) + l] + L, [-5x^ + 5x + L,. (-Sx^ + 8x - 4) + L,. (Sx^ 



8x + 4)-2]+ 6x2 -7x+- l2+ (-10x2 + 1 lx + Lj,(- 12x2 + 12x- 6) _4)l^ + 



73x2 



73x 



■+7I^ 



13x2 13x 13\ 2f.2 7\ / ^ 7 ^ 531 

+ — --j+L2(^6x2-7x+-J+L,(l0x2-9x + 3) + -| 



+- 



X^ X 1 

" 3"^3 ~ 6 



+ L„,l2[2x2-2x+1]+l3 



2x2 2x 1 
"~^y ~ 3 



+ L2 [-x2+X 



-Lj, (-2x2 + 2x - l) + ^2;c2 _ 2;^. _^ _ _ _^ L^^L^ |-_g^2 _^ 5^ _^ ^ (_ 8;(;2 + 8x - 4) 

-L;,(8x2-8x + 4) -2] +l2 [5;c2-5;c + L;,(-8x2 + 8x-4) +Lv(8x2-8x + 4) +2] 

7x2 7^ 7 



- -xj l2 + (3x-2x2) L, + ^ +L2 -xJ - — +7l2 



6 12 



-Ly(2x2-x-l) + 



25 



+ L, [ ( - 1 2x2 + 14x - 7) l2 + (20x2 - 22x + Ly (24x2 _ _^ ^2) 



+8) - ^ + ^ + L„ (-20x2 + 18x - 6) +l2 (-12x2 ^ ^4^ _ 7) + 7r2 - 

13 \ 53 " z' 2jt^ 3 \ z' 3 \ 

+— J - Y + f — + 3x-- jL^.+ (-6x + Lj.(-2x2-2x+l) + - jL2 + Li2(x)(-4x2 



+2x- l)L-t + 

2 /nn„2 



o 3\ o ^ . 55x 2/^46x2 31x 31 

3X-2 L? + 3L,-24x2 + _ + .2/ + _ 



10 



429x 



3L^ 429r -r^^^ . 

— - H h 

2 16 16 



32x2 17x 29, 
■^ + --24)+^.^ 



2x2 5-^ 



+ 0+Si,2(x) (4x2 



, , , n , , , 41x . / 46x2 61x 61 \ 13 
-6x + 3) +Li3(x) (4x2 -2x+ 1) | 24^2 _ _ ^ ^^2 / _ + _____ j + _ 



+ 



34x2 28x 14 \ ^ 2831 



x2 X 1 

T~ 4"^ 8 



2x2 2x 1 
~~T^3 



+L L-* 



4x2 4x 2 
'1~^T~ 3 



X X 1 

y ~3 ^6 



2x2 2x 1\ /2x2 2x 
'1" + T"3 J+^Hl~"T 



+ l2 L, [2x2 _2x + L, (-4x2 _^ 4^ _ 2) + (4;c2 _ 4^ _^ 2) + l] + L„,L,2 [6x2 

10x2 lOx 



-6x + U (-8x2 _^ _ 4^ _^ ^g^2 _ 8^ _^ 4^ _^ 2] + L; 

X 1 



3 ■ + 



^ , 16x2 16x 8 
3/+^H^-- + 3 



+ L 



2 4 



1 



5x 



1 



-4j+y+^H-^ +2+2j+^ l--y+3 



2x2 2x 1 



U+^x- 

3 
4 



3x 



5x^ 



16x2 16x 



■ + L2 



2 J 2 ' y\2 

+ UL4(2x-1)l2+(4x2 



10 



-6;c)L,- ^ + _ + l2(2x- 1) +L, (-4;c2 + 2x + 2) +71^ / — - - + - j - — 

73y2 TO 

;i2;c2 - 14a: + 7) + (-20^;^ + 22x + Lj. (-24^:^ + 24;c - 12) - 8) L, + ^ 



+71^ 



13x2 ^ 13x _ 13 ^ (12,2 _ ^ 7^ ^ (20^2 - I8x + 6) + ^ 



.2t 3 



2x-'L 



+ [-X + Lv (^-2x2 +;c - + 1^ l2 + Li2(x) (-4x2 + 2x - l) L, + (^-8x2 + + 7i2 (2^2 



189x2 



+m--x 



189x 



+ 71" 



5x2 2x 1 



2x2 4^ 



2 \ Z' 1 3jc 

- ^ j + S 1 ,2 (x) (4x2 - 6x + 3) + Li3 (^^^ ^4^2 _ 2;^ + 1 ) + l / 8x2 - ^ + ^^2 ( _2;c2 + 



1)4- 



44x2 









[(4—) 





-4x - 2) - 3) l2 + Li2 (x) ( 8x2 _ 4^ _^ 2) L, + (3 - 6x)l2 - 6Lv + 48x2 - 55x + % 



92x2 



62x 31 



+ 12x-6)+L^ 



20 Lv + 3L; 



429x2 429x 



+ ^ +Li3(x) (-8x2 + 4x-2) +Si,2(x) (-8x2 



4x2 lox 5\ 2/64x2 34^ 29 



92x2 



122x 61 , 
— + y|-13 



68x2 28 \ 283' 

^ r^TJ^^~~l6" 



+ 



1 



3x"-— +— l:+ — 



12 24 



5x2 



+y+L,( 12x^ 



13x 13 



3 



■ + 



17 
12 



-7x2 ^ . 



7x 7 



-5x2 



83x 2 



58x2 34^ 3 35 



l2+ (2x2-4x + 2)L^ + 



/1 16x2 89x 89 



125x 



_ + _J-6jL,-48.' + . ^ 



-F7i^ ( 25x2- 13x-F- 1 



1 
4 

3 3x 

4 ~y 

/ 128x2 



164x 82^ 



2479x2 2479x 



+7^-^ )+Lv( -^ + ^ + f^ )+^' 



12 24 



32 

5x2 2x 1 1 



32 



12 



829x2 191x 191 , ^ / X /o 2 
■W-W + T440)+'^^^(^)('^ 



37x 



20x2 29x 29 \ 37 3 
'^^~~~^) ~ Y^x 



-56x + 28) Li4(x) ( 12x2 _ ^4^ _^ 7) + l2 + 
-hSi,2(x) (10x2- Ilx -hU(-16x2-K56x-28) -FLy (8x2 - 8x-F4)) ^lj^^^) (10x2-k3x 
-hL, (-28x2 + I8x - 9) -h Lj. (I6x2 - 20x -MO) - l) Li2(x) ^ ^22x2 _ l u + l2 

'76x2 52x 26 \\ / 265x2 

3 3~^y 



(-10x2 - 3x -h Lv (-16x2 + 20x - 10) l) L, 7l2 



259x 293 \ 
'~6 12"; 



^3 + Lv(48x2-^ + 7l2 



-25x2 -k30x- 



61 



128x2 140x 
— ^ + — ^ 



11 



70\„ 39\ 2 



mix^ 2S5x , ^ J ^ .N. 247\ 1621 
+ ^+(-2x2 + 2x-l)log(2)-— 



48 



16 



(20) 



\ \ x^ X \\ 1 , 
-2+2-4 +1^^^'" 



2x \\ 19 
T ^3 j ^36 



1 

£ 



X^ X 1 

T ^3 ~ 6 



Ix^ Ix 1 
"~^T~ 3 



x^ X 1 
T^3 ~6 



5x 5x 



2;t:2 



-x + 



1 



+71^ 



40;c2 /4-^2 4^ 2\ 37 

'~9~^~9~^^'[~~T^3) "18 



+ 



649x2 



11x2 _ 

^^^"72/^216/^ " 



589 ) 



108 

x^ X 1 
9"~9 ^18 



2 

649x 
l08" 



+ L, 
'2x2 



5x 5x 



+ 1 



10x2 lOx 



2x 1 
T+3 



"4x2 


Ax 1 










~Y 


-T + 3 



8x2 4 
~9~+ 9"~9 



23x2 23x 23 
l^^Ts"" 36 



, 5 

+4X-- 



+ L 



11x2 



X 1 \ 43 
^6 "12 / ^36 



9 



Ux / 8x2 4\ ^ 



ITx 17\ 19 
T8"^36 / 



4x 2 ,^ 



+ ■ 



2x 1 

T~ 3 



20x2 



22x 



+ L„ 
10' 



73x2 73x 



18 18 
91x2 



, Lv + ■ 
3 / 18 



+ " '-6 

17x2 



91x 2 
7^ + ^ V 18 



+ 



13x , / 2x2 2x 1\ 5\ . ^ . , , / 4x2 



56x2 



■ + 8x + 7i2 



25x2 17x 17 \ 29 , 
^+18 -^'^ 



2 / 2x 



9 V 9 

2 14x 137 \ , , /4x2 4x 2 

+ :rrT +L13W — - — + T I + 



2639x2 2639x 

+ 



V 27 27 216 7 



V 3 



216 

73x2 73x 73 \ 
^ + ^"18]^^ 



216 

3937 



432 



(21) 



2x2 2x 1 
"~^T~ 3 



4x2 4^ 2 
"~^T~3 



10x2 lOx 51 1 , , 



+- 



1 

6 

50x 



+ l2 [2x2-2x+1] +L„ 



5x2 5^ /4^2 4^ 2\ 1 



T ^8x2 



8x 4 

y+3 



19 
"9" 



+ ■ 



131x2 13 Ix 



27 



27 



20x2 



20x 10 



2 

X X 

-y+3 

50x2 



x2 X 1\ 101 

+71" I \ H 

' 3 3 6/ 54 



2x 1 

-T + 3 



14x2 
~9~ + ~9 



14x 7 



"7x2 


7x 


7 ■ 




"4x2 


4x 2" 










+ l2l 








~9~ 


-T + 


18 






-y+3 





2x2 



37x2 37^ / 4^2 4^ 2\ 43' 

1^ iJ^^H r^T~3j^36 



4x2 4^ 



~l~Lm 

14 



9 
76x2 
~9~ 



■ + L, 



76x 



8x2 g^ 
~^T~ 3 J ^9 



+ L 



7x2 



7x , 7 ^ 1 

y + L,(-4x2+4x-2) + - 



20x2 20x 16 \ 9 
~9~^~9~~"9"/ ^2 



+ L, 



40x2 



46x 
~9~ 



199x2 199x 2/^2 1 A 349 



2x2l? / 




' 2x2 


2x 








+ T 











12 



+ 



27 J 

^ 76^2 



4^2 4x 



224^2 272x 7/8^2 16x 



27 



+ ■ 



27 



113^2 ll3x 



■ + 



54 
16x 
~9~" 



54 



17x2 



41x 83 \ , , /4x2 4x 2 
54+T08J+^^^^"H--T + 3 



38 
"9" 



149 
108' 



■ + ■ 



(22) 



X X 1 

2"~2'^4 



+ ■ 



1 



3 

37 



5x^ 



4x2 4^ 2 
5x 



1 



3 

+ L, 



■ 2 

X X 

y ~ 3 

19 

"36 
40x2 



+ ■ 



+ L, 



X 



X 1 

3 + 6 



5x2 



4x2 4^ 





"2x2 


2x 1' 




2 1" 










-X +X-- 




3 


~T^3_ 







40x 



18 



649x2 649x 



108 



108 



20x2 



20x 



8x2 

- + y-3)+L.. 



8x2 4 
"T~y ^3 



4 +n' 



59x2 59^ 59 



20x2 



■+4 





" x2 X 1 ■ 




2x2 


2x 1" 




""9^9 "Is 


+ L2L 


3~ 


+ T-3 





"8x2 


8x 4" 






~9~ 


~y + 9 


+ L2 



55x 



23x2 23x 23 
1^ 18"^36 



4x - 4x + ; 



-l-T T 
5^ 



16x2 16x 8 , 



x2 X 1 \ 43 
+^ '--6 + 12 -36 



16x2 



4x 



16x 8\ 23 
~ + 3 J -9" 



73x2 



4x2 4^ 

-- + T 

55x2 



73x 
18 +T8" 



40x2 44x 20' 



+L 



2 
3 

„2i 3 



4x 

T 



1027x J 



89x2 g9^ g9 
~l8~ + T8"-36 ^' 



26x 



4x2 4^ 2 



10 



-x2L^+ — + LJ — -- + - - L2 + Li2(x) — - — + - L,+ 



40x2 16 



8x2 g^ 4 



1027x2 
54 
787' 
108 

/1 12x2 



+ 



-16x + 7i 
+Li3(x 



50x2 34^ gg^ 



13963x2 13963X . /26x 16 
+ ... +l2| 



8x2 g^ „ / X 

■- + y-3)+Si,2W 



8x2 

-^-3 



648 
8x 4 
3 



648 



3 3 



/ 250x2 202x 631 \ 
27 27~ + 216y 



+ 



13x2 13x 13 \ 



112x2 go^ 



+ ■ 



+ 71^ 



y 

50x2 



4x2 8x 4 
~9~~ 9" + 9 



22x 1 1 \ 26 



10069 
1296 ' 



(23) 



1 



"2x2 



2x 1 

T + 3 



+ L, 



4x2 4^ 2 

~-T + 3 



+l2 [-2x2 + 2x- 1] +L, 
"50x2 50;^- 



-L, 



"5x2 


5x ^ / 


' 8x2 


8x 

+y 


9 




16x2 


16x 8^ 

+ 3y 


1 +L, 


/I6x 


3 


I 3 



10x2 lOx 5 
+ 9~ + ~9~~ 9 



+ ■ 



1 



£ 

8x2 



16x 8' 
- + 3l+- 





2 

m 


r 2 

X X 

J-3 


1" 

+ 6 


8x 




1 " 




3 


18_ 


19" 




131x2 








27 





13 



Ulx 

101 
'54 



AOx^ AOx 20 \ . /5;c^ 



5;c 5\ /40x2 40;c 20' 



7x2 7^ 7 ■ 
"~9~^ 9"~T8 



-l-T T 



4x 2 
'~^T~ 3 



14x2 _ 14x 7 
~9 9~ + 9 



37x2 2,lx 



8x- 



8x 



2x2 2x 1 
'~^T~ 3 



4\ ^ /8x2 



_43 
~36 

+L, 



65x2 



4x2 4^ 

\ h L,, 

9 9 • 



16x2 16x 8 , ^ 



16x2 16x 8 
~i 3~^3 



8x 4 
T + 3 



^ + ^ + L, ( - 8x2 + 8^ _ 4) + L, ( 8x2 - 8x + 4) - 



40x2 40x 32 \ /40x2 40x 32 



+ L, 



2 4x 



76x2 76x 



80x2 



92x 28 \ 451x2 



/ , 2 4x 



45 Ix 



27 



-5x^ + 5x- - ) +Lj, 



80x2 



4 2t 3 /^26x ^ 

+ --x2l;+ — 



68x 16 \ 301 
%x 4\ / 448x2 544x . 

T + 3j^-^+(^-^+" 



4x2 



4x 2 

y+3 



10\ , ^ . , , /8x2 
- L; + U,(x) - 



16x2 32x 16\ 260' 



5809x2 
162 



+- 



5809x 



162 



26x 16 



8x2 §^ 4 



+^{— -ir)+Li3W(-— + --^)+Si,2W(-— + --3 



8x2 4 



4x2 8x 



4\ 487x2 39 Ix 301 \ 

+ - +71^ \ + L,, 

9) V 54 54 108; • 



448x2 352x 
+ ■ 



27 



27 



+ 71" 



16x2 



64x 32 



164 

In 



+ 



16x2 16x 8" 

^ + 9 ) ^3 



5023 
324 ' 



(24) 



4x2 4^ 2 
~9~~ 9"^9 



40x2 40x 20 

+ 



27 27 27 



100x2 lOOx 2/ 4x2 4^ 2\ 50 
^~81 ST^'^ I 9~^9"~9j^8T' 



(25) 



Notice that our results are expressed not only in terms of the classic polylogarithms up to weight 
four, but also in terms of Nielsen polylogarithms 



{n-\)\p\ Jo 
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Conclusions 



In this letter, we have presented the two-loop virtual QCD corrections to the production of heavy- 
quarks in the light quark annihilation channel in the ultra-relativistic limit. Our results form a 
crucial part of the NNLO predictions for heavy-quark production in hadron-hadron collisions. 
However, depending on the mass of the heavy quark (bottom or top) and the kinematics of the 
process under consideration power corrections in the heavy-quark mass may have to be considered 
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as well. Our results have been derived by combining two completely different methods which have 
substantial overlap. This provides direct and highly non-trivial checks of the QCD factorization 
approach of Ref. [13], of the direct calculation of massive Feynman diagrams and, last but not 
least, also on the available masseless results of Ref. [14]. 

In order to yield physical cross sections, our result for (iW ^^^9V[ (^)) still has to be combined 
with the tree-level 2^4, the one-loop 2 — *• 3 as well as the square of the one-loop 2^2 pro- 
cesses. While some of the matrix elements (including the full mass dependence) can be easily 
generated, others became available in the literature only rather recently, see e.g. Refs. [22,38]. The 
combination of all these contributions enables the analytic cancellation of the remaining infrared 
divergences as well as the isolation of the initial state singularities which need to be absorbed into 
parton distribution functions. This is a necessary prerequisite e.g. to the construction of numerical 
programs which provide NNLO QCD estimates of observable scattering cross sections. 

Finally, while we have chosen to work in this letter with squared matrix elements, it should, of 
course, be clear that with the help of Refs. [39,40] analogous results can be derived for the massive 
two-loop amplitude | fAf ) itself. 

A Mathematica file with our results can be obtained by downloading the source from the 
preprint server |http : / / arXiv . org[ The results are also available from the authors upon request. 
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